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The present paper describes a simple guidance strategy for the reentry mission of a 
generic Reusable Launch Vehicle. The guidance computation starts with a nominal feasible 
trajectory, which gets updated during actual reentry based on onboard measurements, using 
both incremental angle of attack and bank angle parameters. An analytical formulation for 
a modified predictor-corrector algorithm is developed, using angle of attack profile and 
bank angle reversal strategy. The guidance solution is arrived at under the heating rate and 
equilibrium glide constraints and the control actuation is implemented through the 
aerodynamic control surfaces. A scheme to improve the overall performance of the 
predictor-corrector algorithm is evolved, which is based on the determination of optimum 
perturbation levels of control variables. In addition, the control vector formulation is 
provided with a weighting option, for improving the overall terminal interface accuracy. 
Numerical results are obtained for the time histories of control variables and performance 
parameters, for a typical nominal trajectory, which show that proposed guidance scheme 
works well for a specified reentry mission. Results are also obtained for the optimum 
perturbation levels for the control parameters, which are used for estimating the sensitivity 
coefficients required in corrector algorithm. Lastly, the results are obtained for the 
trajectory parameters, by considering dispersions in the aerodynamic drag and lift 
coefficients, for establishing its robustness. The results clearly bring out the fact that the 
proposed reentry guidance algorithm is able to achieve the required terminal conditions, 
even in the presence of large perturbations in the aerodynamic coefficients. 


Nomenclature 

ap, 4.5 2pm atm : Accelerations due to Drag and Lift, and Corresponding Measured Values 
C : Heat Capacity 
Cr, Cp : Lift and Drag Coefficients 

c : Path Constraint Function 

a : System Dynamics 

g : Terminal Constraint Function 
h, he : Altitude, The Final Value 
J : Jacobian Matrix 
K, Ky, Ko, Kp, Kp : Various Gains used in Guidance Algorithm 

0,0, Qvnox : Heating Rate, Its Maximum Value 

r : Vector Distance of Vehicle from Centre of Earth 

R, Rw, Rig : Range Solution, Cross Range and Down Range 

ie Ais : Control Vector, Its Initial Value, Its Optimum Calculated Value 
Vi, Ve : Reference Velocity, Relative Velocity 
[w] : Weighting Matrix 

X, Me X; : State Vector, Its Initial Value, Its Final Value 
a, 6, Aa, Ao, o4 : Steering Angles; Control Variables 
P, PsL> Pm : Atmospheric Density, Its Sea Level Value, Its Measured Value 
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I. Introduction 

Reusable Launch Vehicles (RLVs) are being developed as low cost space vehicles, involving the recovery 
of one or more of the stages, after completion of the ascent mission. In view of the fact that the recoverable stage is 
required to de-orbit, reenter the earth’s atmosphere and land at a pre-designated landing site, without violating the 
explicit and implicit constraints, such vehicles require a reliable navigation, guidance and control system which is 
also robust. In addition, it is necessary to make the algorithms, fairly generic as well as re-configurable, in order to 
minimize the cost of operation for a large class of RLVs. It has been established since the developments in space 
shuttle programme’ that atmospheric reentry is the most critical part of the overall return mission and the main task 
of entry guidance algorithm in this phase is to steer the vehicle through the atmosphere, while dissipating the large 
orbital energy, for meeting the mission requirements. 

The flight proven entry guidance algorithms have attempted to do this by fitting the optimal trajectory in a 
phased manner, so as to follow a pre-computed reference trajectory profile’”. However, such guidance algorithms 
are required to be (1) computationally efficient because of the limited onboard computer resources available and (2) 
fairly realistic for ensuring required accuracy at the terminal phase. There have been strategies to address this 
problem, based on analytical expressions for the reference profile, which keep the overall execution time small, but 
suffer from a lack of accuracy and operational flexibility, since an onboard trajectory redesign is restricted to the 
immediate vicinity of the pre-defined reference profile. Current research in this area has focused on autonomous 
guidance, which is based on predicted future paths, as a solution of a Two Point Boundary Value Problem (TPBVP), 
wherein conditions at both starting and terminal points are specified. In theory, there can be infinite paths that satisfy 
the prescribed boundary conditions and it becomes necessary to use an additional constraint to uniquely determine 
the applicable trajectory solution and, in most cases, this is achieved through minimization of a cost function. This 
reference trajectory is also used as a nominal trajectory for implementation of guidance algorithm. 

It should be mentioned here that there are large uncertainties in the vehicle and atmosphere related 
parameters during the re-entry mission so that there is a need to compute the trajectory onboard the RLV in a real 
time manner’. In literature, the standard predictor-corrector algorithm is considered as an appropriate option for 
predicting fairly accurate trajectories, with the conflicting requirements that while computation time should be as 
small as possible, the mathematical models should be as detailed and accurate as possible. A notable study on re- 
entry guidance is by Fuhry* who has proposed the basic methodology of using error sensitivity coefficients for 
correcting the predictions of the trajectory. He has however, used only bank angle magnitude and time of sign 
reversal as controls to meet the terminal constraints while predefined angle of attack is used to meet heating rate 
constraint. Zimmerman et al” have developed an algorithm, which is a combination of predicted corrector algorithm 
where the bank angle profile is estimated to follow constant heating rate boundary with predefined angle of attack 
profile and a profile follower using Linear Quadratic Regulator (LQR). 

In the above studies, the standard form of predictor-corrector algorithm’ is employed, which weighs 
equally, the down range and the cross range error sensitivity coefficients. Further, it is found that in literature, 
evaluation of sensitivity coefficients is based on perturbation of control variables, whose magnitudes are chosen 
arbitrarily, and as such, cannot ensure the correctness of error sensitivity coefficients, particularly in the critical 
phase close to the terminal point. In most of the algorithms reported in literature, only the nominal angle of attack 
profile is used in the predictor algorithm. It is believed that by incorporating some of these features, the performance 
of predictor-corrector algorithm can be substantially enhanced, in terms of improved accuracy, versatility and 
computation time. Therefore, in the present study, a versatile reentry predictor-corrector guidance algorithm has 
been proposed for a typical atmospheric reentry mission, using (1) incremental angle of attack as an additional 
control vector, (2) optimized perturbation of control vectors for improved sensitivity estimation and (3) weights on 
the range related error coefficients to enhance versatility of the algorithm. 

It is believed that these enhancements would allow the use of relatively higher fidelity mathematical 
models of the vehicle and environmental characteristics’, in order to improve overall reentry guidance accuracy, in 
relation to the analytical techniques. In the present study, a weighting option is added to the predictor—corrector 
algorithm formulation, in order to generate the necessary angle of attack and bank angle profiles, to achieve the 
desired targets with the constraints on heat rate and equilibrium glide. The study also proposes the use of 
accelerometers data to improve the accuracy and robustness of the algorithm. Lastly, the study aims to establish the 
robustness of the proposed guidance and control algorithm, to variations in the nominally assumed aerodynamic 
coefficients. 


II. Reentry Guidance Problem Definition 
The objective of any re-entry guidance system is to steer the RLV to a specified Terminal Area Energy 
Management (TAEM) interface, within the prescribed error bands and with adequate robustness, to operate in a 
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dispersed flight environment. The dispersions of significance to any reentry vehicle trajectory include; dispersions 
on (1) vehicle characteristics (mass / aerodynamics), (2) environment characteristics (density, temperature and 
atmospheric winds), (3) initial entry state vector (velocity, flight path angle and heading) and (4) propagation error 
in navigation state. In order to achieve the mission objectives successfully, the vehicle needs to fly the entry corridor 
defined by the thermal, structural, hinge moment constraints and maneuvering capability constraints. Figure 1 below 
represents the typically available atmospheric reentry corridor. 
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Figure 1. Atmospheric reentry corridor in terms of (a) altitude and (b) drag acceleration 


It can be seen from figure | that there is a limited corridor, in respect of a generic RLV, within which the 
required guidance objective is to be met. The overall system dynamics’ is given by, 


X= f(X,TO), At) 
where, X is the state vector, X0_ is the initial condition vector and u(t) is the control vector. The guidance 


objective in the present study is to estimate the control vector history, u(t), so that at final time, ty, the terminal 


state vector x, satisfies the target conditions and path constraint, as given below. 


a(Xp) = 0 c(X(t),u(t)) < 0 (2) 

Of the various constraints that need to be considered during the re-entry mission, the heat rate constraint is 

important in the initial part where the velocities are large. In addition, the equilibrium glide constraint is active 

throughout the entry mission. In the present study, the physical constraints arising out of equilibrium glide path and 

heating rate are incorporated in the following manner. The constraint imposed on stagnation heat rate is as given 
below. 


1/2 
QO. a a . Vn = Orrax 
V Ry Psi V cir 


The methodology adopted is such that the heat rate is clamped to the maximum allowable value when it 
crosses the limit in order to ensure that heat rate value follows the allowable boundary till it falls below the limit. In 
practice, the heat rate constraint can be achieved either by modulating angle of attack or bank angle or both. While, 
the angle of attack modulation has direct control over the heat rate, it can also result in (1) rapid departures from 
predicted trajectory, (2) exceedance of angle of attack limits and (3) possible loss of vehicle trim and stability. 

Therefore, in the present approach, the heat rate constraint is achieved through bank angle modulation, 
based on the required steering computed by the predictor-corrector algorithm. After the computation of 
instantaneous desired angle of attack and bank angle by the predictor-corrector algorithm to meet the target 
conditions, the bank angle is modulated to ensure the heat rate constraint. 


(3) 
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During the closed loop tracking, the coupling between bank angle and angle of attack is expected to be very 
small. Once the heat rate constraint is violated, the vehicle needs to fly at higher altitude, which requires more lift 
force. The heat rate constraint satisfaction is accomplished by computing the incremental bank angle required to fly 
along the specified heat rate boundary and thus the modulating bank angle is computed by the guidance algorithm to 
ensure higher lift. The equation of motion along the altitude is then given by, 

_ Lceosycoso—Dsiny 


(4) 
m 
Assuming small values for 7, the perturbation equation due to Ao is given by, 
- C,S Ao 
po So (5) 
m 
Incremental lift factor (assuming constant C,;, S , m) required to ensure heat rate limit, is taken as, 
qcos AG = K p(Q. ~Qmax) + K p(Qe — Prax) (6) 
Combining above equations with condition that Q.,,, = 0, leads to 
. CS ae 
h-~—Z~LK p(Q. ~ Qmax) + K pO] =0 (7) 


The gains Kp, Kr can be computed by linearizing eq. (7) in altitude with the assumption that time rate of 
change of V, is small. Next, the expression for heat rate can be written as, 


0. = Ci p (8) 
where the density model is taken as exponential. The final form of altitude dynamics is, 
, d : d ‘ 
erly ae) ae h-Qrae) + KpCV qo ae h] = 0 (9) 
m 


which is a standard second order oscillator. In the present study, the @, = 0.05 rad/s and ¢ = 1 are 
assumed as oscillator parameters, for the computation of gains Kp and Ke. Lastly, the incremental bank angle 
required for the heat rate control involves the computation of Q., which is obtained by defining a first order filter 

for heat rate and using both state and its derivative for bank angle computation as shown below. 
Cos( Ac )= Kp(Q-Qmax) + KrQ (10) 

q 

It may be noted that the above formulation is capable of imposing both maximum heating rate and 
maximum total heat load constraints. In the present study, a heat rate constraint of 60 Watts/cm’ is considered to be 
an appropriate constraint magnitude. The equilibrium glide constraint is concerned with the condition that the 
vehicle should not fly above the equilibrium glide boundary. The constraint is defined mathematically through an 


inequality relation, as given below’. 
2 


Peas - cE Hs] cos(v) <0 (11) 
m r 


The implementation of the two constraints discussed above, is done in a sequential manner in which, the 
equilibrium glide constraints are imposed first and the angle of attack and bank angle solutions are generated. Next, 
these solutions are checked if the heat rate constraint also becomes active simultaneously. In case the constraints do 
not become active, the guidance loop is continued, but if the heating rate constraint also becomes active, then bank 
angle is further modulated as shown by equation (10), to satisfy both the constraints explicitly. 


III. Proposed Generic Optimized Guidance Approach 

The proposed guidance algorithm is based on the standard numerical predictor-corrector approach*, whose 
operation is as follows. The predictor integrates the equations of motion using the assumed control profile, including 
the environment and vehicle models’, to estimate the predicted miss (error signals) relative to the target at the final 
time, in terms of longitude and geodetic latitude, defined at the TAEM interface energy level. These error signals are 
used to update the control vector to meet the target conditions, in the following manner. The longitude and geodetic 
latitude are first converted to a target unit vector in Earth Centre Rotating (ECR) frame’. The corresponding 
mathematical representation is as reproduced below. 
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6, = | Re x (12) 


The target unit vector is then given as, 


cos(¢, )cos(/,) 
lee = cos(¢,,) sin(/,, ) (13) 


sin(@,, ) 
Position miss vector is computed in terms of down range (Rpr) & cross range (Rcr) errors as, 


R=([Rpr Rer]* ue) 

It should be mentioned here that the guidance computations are carried out in four different reference 

frames’ namely; Earth Centered Inertial (ECI) frame, Earth Centered Rotating (ECR) frame, Geographic (G) frame 

and the Wind (W) frame. Next, the sensitivity vector of the range position miss vector R_ is defined with respect to 
the control input wu as, 

OR /du ={[OR pp AR |/du} 7 (15) 

Trajectory predictions are performed using perturbed control vectors to estimate the sensitivities of target 

miss to control vector and these sensitivities are used to compute the required control vector to reach the desired 

final target. In the present study, the steering is done through an incremental angle of attack (q) and bank angle (o) . 


In addition, the proposed methodology makes use of the in-flight measurements to increase the accuracy of 
predicted trajectories by compensating for off-nominal flight conditions. The target conditions used in the present 
study are; specified energy level, velocity heading and longitude and latitude at TAEM interface. In the algorithm, 
the predictor is extended up to the required energy level E;. Also, because of the in-built bank angle reversal 
algorithm for steering, heading condition is automatically satisfied. With this, only two constraints of achieving 
defined latitude and longitude are used in the proposed algorithm. Further, the deviations on the angle of attack 
profile are constrained, in order to satisfy vehicle constraints, whereas the bank angle profile is used as the main 
control parameter. Angle of attack and bank angle profile boundaries, considered in present study, are given in 
below in figures 2 (a) & (b). 
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Figure 2. Proposed Profile for (a) angle of attack, @ and (b) bank angle, o 


Here, the mean value of & is considered as nominal angle profile and variation over nominal profile (Aa ) 
is considered as one of the control parameters. Generally, the sensitivity to bank angle is smaller at higher altitude 


than lower altitude. Hence, the bank angle history is selected with a desired bank angle, Oy at a predefined velocity, 
V; and linear variation of the bank angle with respect to relative velocity as given in Figure 2(b). The magnitude of 


bank angle, Og at velocity, Vj; is assumed as another control variable. Bank angle sign reversal strategy is used to 
drive the velocity heading error to an allowable boundary. 
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IV. Predictor-Corrector Algorithm 
Figure 3 given below shows the schematic of the predictor-corrector algorithm that is proposed in the 
present study. The main components of predictor algorithm are; vehicle and environment models, in-flight 
estimation and equations of motion. 


Nominal Trajectory Prediction 
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Target Miss Sensitivity 
Control Vector Correction 
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Figure 3. Schematic of Proposed Predictor-Corrector Algorithm 


The translational equations of motion are expressed in Earth Centered Inertial frame (Z-axis is takes as 
passing through north pole, X-axis along Greenwich meridian at the time of entry and Y axis completes the right 
handed triad) as given below’. 


r=V ; V = g +4,+4, (16) 
where, 7 V are vehicle inertial position and velocity vectors, Z is gravitational acceleration and 4p, Gy are 


acceleration due to drag and lift respectively. The environmental effects considered are; variation of atmospheric 
properties with altitude, Earth’s oblateness effect on gravity vector, effect of atmospheric rotation with Earth and 
movement of target point due to Earth’s rotation, while the atmospheric density, temperature and wind velocity are 
computed from the stored profile. 

The Indian Standard Atmosphere (ISA) is used for the predictor, with the assumption of wind data update 
available before initiating reentry. The vehicle characteristics considered are; vehicle mass, aerodynamic 
coefficients with respect to different flight parameters and control history. The nominal aerodynamic lift and drag 
force coefficients for different aerodynamic flow regimes are used. These coefficients are functions of altitude, 
viscous interaction parameter and Mach number depending on the aerodynamic regime. For computational ease, 
these tabular data are stored as curve fits and the nominal aerodynamic coefficients are computed using these curve 
fits. Based on the in-flight estimates, the aerodynamic accelerations are updated as given below: 


dy = (1/2)K, Pp V_ [Cp/m| Sip 3 4 = 0/2) Kes py (L! Dap (17) 
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where, K p> K(z/p) are estimator coefficients, V, is relative velocity, § is reference area, 9 is nominal 


density, mis vehicle mass and Cp, L/D are nominal aerodynamic parameters. Uncertainties on vehicle and 
environment parameters are compensated in the predictor by applying suitable scale factor, which is estimated, 
based on actual acceleration at any point in trajectory. Using the accelerometer data, the measured drag (ap, ) and 


lift (az,,) accelerations are computed. The primary dispersion affecting the aerodynamic acceleration is the 


atmospheric density and hence, for implementation simplicity, the dispersions in the aerodynamic acceleration due 
to atmospheric uncertainties are lumped into a density scale factor as, 


Kp =Pm!P: p, =(2a,,/V,)|mi(C,S )| (18) 


The nominal vehicle characteristics and nominal density are determined using the predictor models for the 
vehicle state at the time of measurement. Since nominal ballistic coefficient is assumed for deriving the measured 
density and measured acceleration is based on the actual ballistic coefficient, uncertainties in the ballistic 
coefficients are also reflected in the measured density. The uncertainties on lift coefficient can be lumped into a 
scale factor given by, 


Kp) =L/D)m\/[(L/D)) ; (L/D) m= 4tm!/4Dm (19) 
Measured accelerations are noisy signals and exhibit short-term variation due to short-lived local 


atmospheric dispersions. Filtering of scale factors is therefore necessary and implemented using first order filter as 
given below. 


Kp =(1—K1)K p + Ki(Pm/ p) (20) 
K(t/p) = (-K2)K(z/p) + Ka[(L/D)m KL/D) (21) 
where, Aj, K,are filter gains. The proposed re-entry phase guidance algorithm controls the two final state 
constraints viz.; down range error and cross range error, using the following control vector, 
— T 
u=(Aa@ og] (22) 
In earlier studies*, a Taylor series expansion of the target range error, about the reference control solution, 
u » 1s proposed, for generating control vector, as follows. 


RW) = Ro) + Jot - 9) + (23) 
iw =i9 — J~' (up) R(t) (24) 
Sos ci 
where, the Jacobian matrix is defined as / (uo) = Oat (25) 
Uo 


The above formulation for control vector is limited by the fact all the error sensitivities are weighed 
equally. In order to make the re-entry guidance algorithm more generic as well as re-configurable for different kinds 
of re-entry vehicles and missions, the present study proposes to include a weighting matrix for the range error 
sensitivities, as follows. 

* 


— = ey an 
uo =u,-—J (u,) [Ww] Ry) (26) 
Here, [w] is a matrix of weights that can be selected by user a priori that are suitable for a specific re-entry 


vehicle and mission. When all elements of [w] are unity, it reduces to the formulation, as described by equation 
(24). It may be noted here that the partial derivatives required for the generation of Jacobian matrix are not available 
analytically and in the proposed corrector algorithm, similar to the methodology adopted in other references*”, these 
are computed numerically onboard, using the results of three trajectory predictions performed with three different 


control vectors viz., nominal control vector, control vector with perturbed A@ and perturbed Cd . However, the 


success of numerical predictor-corrector algorithm lies in establishing the correct sensitivities of the state errors, 
relative to the control variable. Further, in order to implement such algorithms on-board using limited computational 
power, it is necessary to obtain correct sensitivities within least possible time. Such an approach has the potential to 
avoid the need to generate hundreds of trajectories to establish the necessary functional relations and to store them 
on-board for use during the mission. 
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It should be noted here that the correct selection of perturbations of control variables for generation of 
numerical derivatives is crucial for the success of predictor-corrector method. It is well known that if the 
perturbation level is too small, the derivative is meaningless because of loss of accuracy due to round-off errors. On 
the other hand, if the perturbation level is too large then the process will give only an averaged, and approximate, 
value of the required derivative. 

Therefore, in the present study, a methodology is evolved for obtaining an appropriately correct value of 
the perturbations in the control vector, as a function of the vehicle relative velocity, which can be used for obtaining 
the correct error sensitivities on-board the vehicle and in real time. The methodology involves carrying out a large 
number of off-line simulations on ground and determining the perturbation size that minimizes the two errors of 
interest. These perturbation sizes are then stored on-board, in the form of look up tables, as a function of the vehicle 
relative velocity. It is believed that such a methodology provides a solution that can be adapted to a wide variety of 
RLV configurations and missions. It should be noted further that, (1) Optimum perturbation level at each instant of 
flight is different, (2) Perturbation levels have unique relationship that is valid for a range of flight environment, (3) 
Perturbation level and its variation pattern are different for different control variables, (4) Different levels are needed 
for entry and terminal phases, for the same control parameter and (5) Sensitivity coefficients depend on flight 
condition, for the same perturbation level. 

It is to be kept in mind that, even with better estimates for the sensitivities, the corrector algorithm provides 


only an approximate solution and hence the control vector solution 7” is achieved through an iterative process. The 
sign of bank angle is chosen to reduce the heading miss, which is defined as the angle between current plane of 
motion and the plane formed by current position vector and target vector and is based on the sign and magnitude of 
heading miss. Towards the end of trajectory, a small error in the target miss can demand large attitude commands 
and it is desirable to minimize such maneuvering demands. Hence, a terminal guidance law’ is introduced which 
minimizes the target miss and avoids instability. In this case, it is assumed that only the bank angle o is available in 
the terminal phase. The guidance problem then involves solving for the bank angle increment Ao over the profile 
predicted by the entry guidance that results in the minimum achievable target miss. The minimum miss distance, for 
the linear approximation, lies along the normal to the line from the target as shown by terminal condition vector 
diagram, given in Figure 4 below. 


A ; “\ Vehicle Position 


AR min, End conditions with 6. 


/ oe Target i’on 
la 4 Ri (ota) 
4 End conditions with (Go+8o) 


Figure 4. Terminal Conditions 


The approximate bank angle corrections required to achieve minimum miss is computed as 


Ao = [ARmin [aR] /|oo|| ' SIGN (AR- ARmin) ; where, AR=R—Ro_ 07 


V. Results for Nominal Reentry Trajectory 
The proposed guidance algorithm is coded as a separate module and integrated with a simulation tool®. The 
guidance functions are called at every cycle and ideal navigation is assumed. The state vector and accelerations of 
the simulator are passed on to the guidance system as navigation information and navigation is assumed to be error- 
free. Similarly, it is assumed that control commands are generated instantaneously and the steering angle computed 
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2. 


HEAT RATE (W/cm ) 


by the guidance system is fed back to the trajectory simulator. Typical data available in literature is used for the 
vehicle characteristics’, whereas standard models are used for simulating the environment. The nominal reentry state 
vector and the target conditions used for the present study is given in appendix. It is to be noted that these 
parameters are not related to any RLV mission and are selected only to show the performance of the developed 
algorithm. Figures 5 to 8 below provide the various results for a nominal trajectory, with and without imposing the 
heat rate constraint. It can be seen from these figures that nominal trajectory nearly satisfies the heating rate 
constraint, excepting for a small part of the trajectory. It can also be seen that when constraint is active, the heating 
rate in that region is clamped to the boundary value, until the heating rate comes down. The associated change in 
altitude can be seen in the corresponding altitude profile. The results in respect of control variable time histories 
show that while angle of attack is only marginally influenced by the presence of heating rate constraints, the bank 
angle is drastically modified in the region where the heating rate constraint becomes active. With regard to the down 
range and cross range solutions, it can be seen that while down range profile remains unaffected by the heating rate 
constraint, the cross range profile shows a marginal departure in the region where the heating constraint becomes 
active. It is further seen that in both the cases the terminal range conditions are achieved with the same level of 
accuracy and that the control variables go to zero at the terminal point. 
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Figures 9 to 12 present the results for steering angles, altitude, range parameters and heating rate time 
histories, when both the heat rate and equilibrium glide constraints are imposed. It is seen from heat rate profile that 
combined effect of heating rate constraint and the equilibrium glide constraint is almost same as the heating rate 
constraint alone, which indicates that heating rate constraint supercedes the equilibrium glide constraint. However, 
the altitude profile shows that when equilibrium glide constraint is active, vehicle initially flies at lower altitudes and 
flies higher altitude when the heat rate constraint is active. Lastly, the profiles for down range and cross range show 
that while down range profile is practically unaffected whether or not the constraint is active, the cross range 
magnitudes are lower when the equilibrium glide constraint is made active. This is so because the equilibrium glide 
constraint indirectly drives the out of plane forces to zero or less, so that the range magnitudes are lower. Thus, it 
can be seen from the results presented for the nominal trajectory that while heating rate constraint becomes active 
for only a small duration of the re-entry flight, whenever it becomes active, it supercedes the equilibrium glide 
constraint. Further, influence of these constraints is more on bank angle profile and cross range and marginal on 
altitude, while other parameters remain practically unaffected by applying both these constraint. 
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VI. Numerical Results for Optimized Sensitivity Coefficients 


Figure 13 below shows the variation of optimum control perturbation levels as a function of the relative 
velocity that minimizes the error in sensitivity estimation. It can be seen from the figure that towards the end of the 
trajectory, it is necessary to use much larger perturbation levels, in order to ensure sufficiently accurate 
determination of derivatives. Figure 14 below shows the variation of the down range and cross range sensitivities 
with respect to the bank angle, as a function of the relative velocity and it can be clearly seen that as velocity reduces 
towards the end of the reentry phase, the down range variation is practically insensitive while cross range variation 
shows a small but finite sensitivity to bank angle parameter. This is so because lower dynamic pressure requires a 


higher bank angle to achieve the same control force. 
300 
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Figure 13. Optimum perturbation levels 


Figure 15 shown alongside, gives the 
variation of the down range and cross range 
sensitivities with respect to the angle of attack, as a 
function of the relative velocity and it can be noted 
that as velocity reduces towards the end of the reentry 
phase, the down range sensitivity approaches a 
constant value while cross range variation becomes 
practically insensitive to the incremental angle of 
attack parameter. Lastly, figures 16 given below, 
show the optimum perturbation level in bank angle 
and the corresponding total range sensitivity at the 
terminal phase and it is found that optimum 
perturbation levels increase linearly with decreasing 
velocity, while range sensitivity reduces and settles to 
a small finite value as terminal phase is reached. 
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Figure 14. Downrange & crossrange sensitivities 
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Figure 15. Downrange & crossrange sensitivities 
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Figure 16. Optimum perturbation level of bank angle and terminal range sensitivity 


VII. Robustness of Guidance Algorithm to Aerodynamic Parameter Variation 

Next, solution is obtained for the control variable time histories as well as time histories of the trajectory 
range parameters i.e. downrange and crossrange, for +10% dispersions on both drag coefficient Cp and lift 
coefficient Cy. Figures 17-22 show these results and it can be seen that while angle of attack time history shows 
some effect of the dispersions on both C, and Cp, the bank angle history is more or less unaffected by these 
variations. It is also seen that in all the cases (figures 17 & 18), the control command approaches zero towards the 
terminal phase of the reentry guidance. The results presented in figures 19 & 20 show that, dispersions of +10% in 
both lift and drag coefficients, have practically no effect on the actual realization of the specified downrange and 
crossrange. In addition, the total time for the reentry mission remains the same. Lastly, figures 21 & 22 show the 
effect of dispersion in aerodynamic coefficients on the altitude time history and it seen that in case of dispersion on 
drag coefficient, only the terminal part of the time history is affected and also the required altitude is reached later if 
the drag coefficient is lower than nominal and vice versa. It is further found that in case of dispersion on the lift 
coefficient, there is also a noticeable effect on the complete altitude time history, in addition to the above effect, as 
seen for the drag coefficient dispersion. These results presented in this study have clearly brought out the fact that 
the proposed guidance scheme, using bank angle and angle of attack as control parameters, is fairly robust and can 
be employed for a large class of RLVs, having significantly different aerodynamic coefficients. 
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Figure 19. Ranges with dispersions in Cp 
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Figure 21. Altitudes with dispersions in Cp 
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VIII. Conclusions 
A generic predictor—corrector algorithm for reentry guidance of RLV missions is presented in this paper. 
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The complex reentry guidance problem is formulated as Two Point Boundary Value Problem and constraints on 
heating rate and equilibrium glide are imposed on the solution. The guidance algorithm uses incremental angle of 
attack, in addition to the bank angle, to achieve the control objective. The control vector estimation, in proposed 
predictor-corrector algorithm, is also provided with a weighting matrix, to allow simulation of a wide variety of 
RLV configurations and missions. The accuracy under dispersed flight environment is enhanced with the use of 
appropriate estimated parameters in the predictor algorithm. Sample results are obtained for verifying the algorithm. 
The sensitivity analysis establishes a mechanism for estimating onboard, the optimum perturbations for specific 
functional relation, instead of storing a large number of dispersed trajectory solutions. Lastly, a robustness check is 
carried out with respect to dispersions in the aerodynamic lift and drag coefficients and it is found that algorithm is 
fairly robust against significantly large dispersions. Therefore, the proposed guidance scheme is considered to be 
fairly generic and applicable to a variety of entry missions and vehicle configurations. 
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Appendix: Reentry Interface and Target Parameters 


The nominal entry state vector considered for this study is given below: 


Altitude :120km Latitude : 0.0° 
Longitude :0.0° Inertial velocity : 8086.0 m/s 
Flight path angle :-3.2° Flight azimuth angle : 100.0° 
The target conditions at the Energy = 1.2344x10° m’/s’ for terminal phase are: 
Geodetic latitude :-17.98° Longitude : 45, 35° 
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